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In this paper, we study complex Wishart processes or the so-called Laguerre processes (Xt)t>o- 
We are interested in the behaviour of the eigenvalue process; we derive some useful stochas- 
tic differential equations and compute both the infinitesimal generator and the semi-group. 
We also give absolute-continuity relations between different indices. Finally, we compute the 
density function of the so-called generalized Hartman- Watson law as well as the law of 
To := inf {t, det(Xt) = 0} when the size of the matrix is 2. 
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1. Introduction 

The Real Wishart process is a symmetric matrix-valued process which was introduced 
by Bru [2] as follows. Let B t = {Bij(t))ij be an n x m Brownian matrix and define 
X t =BjB t . The process (X t )t>o satisfies the following stochastic differential equation 
(SDE): 

dX t = Bj dB t + dBjB t + nl m At 

= ^X~tdN t + dN?^X~ t + nl m dt, X a = B^B , 

where I m denotes the unit matrix, the superscript T stands for the transpose, \pC t is the 
matrix square root of the positive matrix Xt and (Nt)t>o is an m x m Brownian matrix. 
This process is called the Wishart process of dimension n, of size m, starting from Xq, and 
is denoted by W(n,ni,Xo). Then W(6, m, Xq), where 5 runs over the Gindikin ensemble 
{1, . . . , m — 1} U }m — 1, oo[, is defined as the unique solution of the latter SDE with S 
instead of n. Thus, it can be viewed as an extension of the squared Bessel process to 
higher dimensions. In this way, Donati-Martin et al. [6] tried to extend some well-known 
properties of the squared Bessel processes to the matrix case and derived expressions such 
as the Laplace transform and the tail distribution of some random variables, in which 
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many multivariate special functions with symmetric matrix arguments appear, such as 
gamma, modified Bessel and hypergeometric functions Muirhead [20]. However, the latter 
case is quite complicated to deal with and, to our knowledge, there are no more precise 
results on the law of these variables. Nevertheless, in the complex case, hypergeometric 
functions with Hermitian matrix argument can be expressed as a determinant of a matrix 
whose entries are one-dimensional hypergeometric functions. In fact, Gross and Richards 
[10] established the following result: 

_ det(x™~ J p^ q (ai - j + 1, . . . , a p - j + 1, . . . , b q - j + l;xj)) 

V(X) 

where X is an m x m Hermitian matrix, (xi) are its eigenvalues, p ^ q denotes the stan- 
dard hypergeometric function with scalar argument, V(X) = Y\i<j( x i ~ x j) i s the Van- 
dermonde determinant, and p F q is the hypergeometric function with Hermitian matrix 
argument defined by 

V (n n h h ywrV {ai)r ■ ■ ■ (a P )r C T (X) 



where r = (ri, . . . ,r m ) is a partition of length at most m and of weight k (i.e. ti > T2 > 
• • • > T m , X)i T i = k)^ ( a )r is the generalized Pochammer symbol and C T is the so-called 
zonal polynomial. We refer to Macdonald [19] for further details and Lassalle [16, 17] 
for analogous expressions for multivariate orthogonal polynomials. The detcrminantal 
representation above is due to the fact that the zonal polynomial is identified with the 
(normalized) Schur functions defined by: 



det(*r J ) 



-(xi,.. 



Consequently, one can use integral representations as well as other properties of stan- 
dard hypergeometric functions to obtain, at least when m = 2, some results hitherto 
unknown in the Wishart case. The rest of this paper is organized as follows. In Section 
2 we introduce the Laguerre process of integer dimension. In Section 3 we study the 
behaviour of the eigenvalue process. Then, in Section 4, we define the Laguerre process 
of positive real dimension. Section 5 is devoted to the absolute-continuity relations, from 
which we deduce the Laplace transform of the so-called generalized Hartman-Watson 
law as well as the tail distribution of To, the first hitting time of 0. In Section 6 we focus 
on the case m = 2 for which we invert this Laplace transform. In Section 7, wc compute 
the density of So := 1/(2Tq). Section 8 presents our conclusions. Some special functions 
of interest are discussed in the Appendix. 
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2. Laguerre process of integer index 

Let B be a n x m complex Brownian matrix starting from Bq; that is, B = (-By ), where 
the entries Bij are independent complex Brownian motions, so we can write B = B 1 + iB 2 , 
in which B , B 2 are two independent real Brownian matrices. We are interested in the 
matrix-valued process X t := B\B t which satisfies the following SDE: 

dX t = dB* B t + B* dB t + 2nl m dt. ( 1 ) 



Definition 1. (X t )t>o is called the Laguerre process of size m, of dimension n and 
starting from X Q = BqB , and will be denoted by L(n, m, X ). 

Remark 1. For m = 1, (X t )t>o is a squared Bessel process of dimension 2n, denoted by 
BESQ(2n,X ). 

Remark 2. Set X t = (Xij(t))ij. One can easily check that 

dX u (t) = 2 V / X u (t)d~f l {t) + 2ndt, l<i<m, 
where (p(i)i<i< m are independent Brownian motions, thus Xa is a BESQ(2n, (Xq)u). 

Remark 3. The equation above implies that 

d(tr(X t )) = 2 Vtr(X t ) d/3 t + 2nmdt, (2) 

where (3 is a Brownian motion. Consequently, (tr(X t )) 4 > is a BESQ(2nm, tr(X )) of 
dimension 2nm starting from tv(Xo). One can also deduce from (1) that, for every 
i,j,k,le {l,...,m}, 

(dX^,dX kl ) t = 2(X u S kj + X kj S u ) dt, 

which differs from equation (1-1-5) derived by Bru [2] since, for a complex Brownian 
motion 7, wc have d(7,7)t = and (7,7)4 = 2t. 

Let H m ,H+ be respectively the space ofmxni Hermitian matrices and the space of 
m x m positive definite Hermitian matrices. On the space of Hermitian matrix-argument 
functions, we define the matrix-valued differential operators 

d _( d \ d _ / d \ d f d , 8 \ 

dx' \i).r, i ) ji ; dy' {tin,!. J ; ,/ dz' \dx jk dy jk J jk 

We also define 

fd_V x ^ d 2 f^_±\ \ - 92 

\dz) ifj ' dz lk dz kj ' \dx dyj tj ' ^ dx ik dy k3 ' 
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Proposition 1 (Infinitesimal generator). Suppose we have functions f satisfying 

df df df df 



for all i,j. 



dyij dyji 

Then the infinitesimal generator of a Laguerre process L(n,m,x) is given by 

d 



££ = 2ntr( ft( — 

az 



az 



(3) 



where d/dz is the operator defined above. 

Remark 4- Using the fact that x T =x, y T 
any two matrices A and B, we can see that 

( AA^l - r (—— \ - ( — — 
\ dy dx J \ dx dy ) \ dx dy 



and tr(AB) = tx{BA) = tr{B T A T ) for 



dz 



dx dy 



3. Eigenvalues of Laguerre process 

In this section, we will suppose that n> m. The following result is due to Konig and 
O'Connell [15], Katori and Tanemura [14] and Bru [1] in the real case: 

Theorem 1. Let Ai(t), . . . , X m (t) denote the eigenvalues of X t . Suppose that at time t = 
all the eigenvalues are distinct. Then the eigenvalue process (Ai(i), . . . , X m (t)) satisfies 
the following stochastic differential system: 



E 

k^i 



X t (t)+X k (t) 



Xi(t)-X k (t) 



(It. 



1 < i < m, t < r, 



where the (/3j)i<i< TO ar £ independent Brownian motions and r is defined by r := 
inf{t, Xi(t) = Xj(t) for some (i, j)}. 

Remark 5. With the help of the SDE satisfied by the eigenvalues, we can compute the 
ones satisfied by both processes (tr(X)) and (det(X)): the first has already been done. 
For the second, we find that for t < Tq := inf{t, det(X{) = 0} and for r G R, 

d(dct(X t )) = 2det(Xt)y / tr(X^ 1 )d^ f +2(n-m+ 1) dct(X t ) tr(X t -1 ) di, 
d(log(det(X t )) = 2Y / tr(X t - 1 )dzy t + 2(n - m) tr^" 1 ) dt, 

d(dct(X t ) r ) = 2r(dct(X t )) r ^tr(Xr 1 )dv t + 2r(n-m + r)(dct(X t )) r tr(Xf 1 )dt, 
so we can see that for n = m, log(det(X)) is a local martingale and so is (det(X)) m ~ n . 
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Lemma 1. Take X a e H+. Then for n>m, X t £ H+. 

Proof. In fact, this result is a direct consequence of the fact that, for n = m, logdet(X) 
is a local martingale, and so is (det(X)) m ~ n . Hence, for n > m, these two continuous 
processes tend to infinity when t — ► To, which is possible only if To = oo, because every 
continuous local martingale is a time-changed Brownian motion. □ 

Corollary 1. If Ai(0) > • • ■ > A m (0), then the process U defined by 
is a local martingale. 

Proof. We could follow the proof given by Bru [1] or make straightforward computations 
using the derivatives of the Vandcrmonde function. But we prefer to use a result from 
Konig, W. and O'Connell [15]: for n>m, the eigenvalue process is the ^-transform (in 
the Doob sense) of the process obtained from m independent BESQ(2(n — m+ 1)). Thus, 
if G and G denote respectively the infinitesimal generators of these two processes, then, 
G(h) = and, for all C 2 function /, 

G(f) = ^G(Vf) G{U) = ±G(1) = 0. D 

Corollary 2. //, at time t = 0, the eigenvalues of X are distinct, then they will never 
collide; that is, r = oo almost surely. 

Proof. This result follows from the fact that the continuous process U tends to infinity 
when t — ► t, which is possible only if r = oo almost surely (We use the same argument 
as before.) □ 

The proof of the following result is similar to the one derived by Bru [2] in the real 
case: 



Proposition 2 (Additivity property). If (X t )t>o and (Y t )t>o are two independent 
Laguerre processes L(n,m,Xo) and L(p,m,Yo) respectively, then the process (X t + Y t )t>o 
is a Laguerre process L(n +p, m, Xq + Yq). 

We now introduce the Laguerre processes of non-integer dimensions 5. 



4. Laguerre processes with non-integer dimensions 

Let X be a Laguerre process L(n,m, Xq) with n > m. If Xq £ H+, and if y/X~t stands 
for the symmetric matrix square root of X t , it is easy to show that the matrix O defined 
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by O t := \PQ B*, where X t = B*B t , satisfies 0*0 = OO* = I m . Thus, 

d 7t = O t dB t = y/xT 1 Bt dB t 

is an m x to complex Brownian matrix. Replacing this expression in (1), one obtains 

dX t = ^/x~ t A lt + dtfy/X~ t + 2nl m dt. 

Theorem 2. If (B t ) is an to x m complex Brownian matrix, then for every Xq G iT+ 
and for all 6 >m, the SDE 

(4) 



dX t = sJX t dB t + dB* y/X t + 26 I m dt 

has a unique strong solution in iJ+ . Furthermore, if the eigenvalues are distinct at time 
t = 0, then they satisfy the stochastic differential system 



dXi(t) = 2^/Xjt)df3i(t) + 2 



E 



X l (t) + X k (t) 



\i(t)-\ k {t) 



(If, 



1 < i < m, 



where the (f3i] 



Ki<m 



are independent Brownian motions. 



Proof. The proof of the second part of the theorem is the same as before, with S instead 
of n. All that remains is to prove the first part. Note first that (det(JQ)), (logdet(Xt)) 
and (det(X() r ) satisfy the same SDE with 8 instead of n. Hence, arguing as before, we 
can see that To = oo almost surely. On the other hand, the map a i— > a 1 ! 2 is analytic in 
if+ (see Rogers and Williams [22], page 134), so the SDE has a unique strong solution 
for all t > 0. □ 

Definition 2. Such a process is called the Laguerre process of dimension 5, size to and 
initial state X . It will be denoted by L(5,m, X ). 

Remark 6. Any process (Xt)t>o solution of (4) is a diffusion whose infinitesimal gen- 
erator is given by 



S£ = 26 tr 3? 



tr x$t 



tr yxs 



Remark 7. A simple computation shows that 

d(Xij , X H ) t = 2{X a {t)6 k3 + X kj (t)6 a ) dt, for alH,j, k, I e {1, . 



'}• 



We now focus on both existence and uniqueness when 6 > m — 1 and Xq € - see 
Bru [2] for the real case. 

If A is a Hcrmitian matrix, let X + be the Hermitian matrix max(A, 0). If we denote 
by (Aj) the eigenvalues of X, then (Xf = max(Ai, 0)) are those of X + . 
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Theorem 3. For all 5 e R + and X Q = x e H m , the SDE 

dX t = <Jx+dB t + dB* <Jx+ + 2SI m dt (5) 

has a solution in H m . 

Proof. The mapping a i— > \fa+ is continuous on H m . Hence, X exists up to its explosion 
time (Ikeda and Watanabe [11], page Theorem 2.3). Furthermore, from 

2 



\\VX+\\' + \\6I\\ 2 <5 Z + \\X\\ 2 < C(l + \\X\\ 2 ), 

we can deduce that this explosion time is infinite almost surely (Ikeda and Watanabe 
[11], page Theorem 2.4). □ 

Proposition 3. J/Ai(0) > • • • > A m (0) > 0, then, for all t < S := inf{i, Ai = Xj for some 
the eigenvalues of X + satisfy the following differential system: 



dA,-(t) = 2 ] /xUT)dMt) + 2 (* + E %W-\l(t) ) dt ' 1 " l - m - 

Proof. This differential system can be shown in the same way as in Theorem 1, using 

(dX-«,dX«) t = 2(X+(t)5 kj + X+(t)S u )dt, for all i,j, k,l e {1, . . . ,m}. □ 

Proposition 4. J/Ai(0) > •■ ■ > A m (0) > 0, then, for all 5>m — l,t>0, X m (t) > 0. 

Proof. First, we note that S = oo almost surely. Indeed, one can easily show that the 
process U defined by 

U(\ 1 (t),...,X m (t)) = 



is a local martingale. For the proof, we proceed along the lines of Bru [1]. □ 

Theorem 4. 7/Ai(0) > • • • > A m (0) > 0, then, for all 8 > m— 1, (4) has a unique solution 
in in the sense of probability law. 

Proof. By Proposition 4, the solution of the SDE (5) remains positive for all t > 0, thus 
it is a solution of (4). □ 

Theorem 5. Let be the space of positive Hermitian matrices. Then, whenever the 
SDE (4) has a solution in H^~ n , for fixed t, its distribution is given by its Laplace transform 

E Xo (cxp(-trwXt)) = (dct(/ m + 2tu)y S cxp(-tiiX Q (I m +2tu)~ 1 u)), (6) 
for all u in . 
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Proof. For s £ H+, let g(t, s) = A~ s exp(-V(i, s)), where 

A t =dct(I m + 2ut), W t ^(I m + 2ut)- 1 u, V(t,s)=tT{sW t ). 
First, note that W G H m . To proceed, we need the following lemma: 

Lemma 2. The function g satisfies the heat equation: dg/dt — J£g, where Jz? is the 
infinitesimal generator of X . 

Proof. If we write s = x + iy, then, using the fact that x is symmetric, y is skew- 
symmetric and W is Hermitian, we can see that tr(sWi) = tr(xM + iyN), where 

W + W Ar W-W 

M= , N = . 

2 2 

Observing that M T = M and N T = -N, we can deduce that g satisfies the conditions 
of Proposition 1. Furthermore, 

^ = -g(28 tr(^)-2tr( S VF t 2 )) 
at 

tr ( x (S " SO) = gt ' (x{M ' + N ' }} = gtI ' lW '' ) - 

Finally, noting that tr(M) = tr(W), wc obtain the equality. 

Now, we consider the process (Z(t, X t )) defined by Z(t, X t ) = gif \ — t, X t ) for all t < t\ 
for fixed t\. From the lemma, we deduce that Z is a bounded local martingale and thus 
is a martingale. So, the result follows from a simple application of the optional stopping 
theorem. □ 

Corollary 3. Let (X t )t>o be a Laguerre process L(S,m,x) where x € H^. For 5 > m— 1, 
its semi-group is given by the following density: 

Pf (*, V) - exp (-1 tr(x + y)) (dct ^ (s; l {y>0} 

with respect to Lebesgue measure dy = Y[ p < q ^Vpq llp< 9 ^Vpq where y = y 1 + iy 2 and qFi 
is a hypergeometric function with matrix argument ( Chikuze [4J; Gross and Richard [10]). 

Proof of Theorem 5. In fact, this result can be easily deduced from the case where 5 = 
n is integer, since, in this case, X t is a non-central complex Wishart variable W(n, 2tl m , x) 
James [12] with density given by 

/«(», V) = {2t)m ! Tm(n) «P + V)) (** VT- (n; f ) l {y>0} 
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with respect to dy. Hence, taking 6 instead of n and denoting by Wt this new variable 
(starting from x), we can see that, using |y| to denote det(y), 

E ( e -truW t ) _ 1 e xo(- — 

x[ > (2t)™ s T m (S) P \ 2t 

x j cxp (-1 tr ((7 + 2 «t)y)) |y| 5 -'" Fx (<J; dy 

2^ m,5 |a;h ,5 / trx s 
■ exp 1 



r m {5) "V a 

x / cxp(-2ttY{x~ 1 / 2 {I + 2ut)x- 1 / 2 z))\z\ s " m Q F 1 (S;z)dz 



z>0 



| J + 2ut \~ s cxp ^- ^ tr(a;(/ + 2ut)- 1 (I + 2ut - 
\I + 2ut \~ 5 exp ( - tr(x(I + 2ut)~ 1 u) J , 



which is equal to (6). □ 

Remark 8. In the last proof, we used the change of variables z = x x / 2 yx 1 / 2 which gives 
dz = |a;| m dy. For the second integral, see Faraut and Koranyi [8], Proposition XV. 1.3. 



Remark 9. The expression for the semi-group extends continuously to the degenerate 
case, namely: 

Pli^v) = {2t) J Tm{§) ^(-^f 1 ) (dety)'-l {l>> o } , 
where m denotes the null matrix. 



Corollary 4. For 5>m—l, the semi-group of eigenvalue process is given by 

where x = (x 1} . . .,x m ),y = (yi, . . . ,y m ) so that x x > ■ ■ ■ > x m > 0, yi > ■ ■ ■ > y m > 0, 
5 = m + v such that v > — 1 and l v denotes the modified Bessel function Lebedev [18]. 



Proof. The expression for the semi-group can be computed using Karlin and MacGregor 
[13] formula since, for 5 > m — 1, the eigenvalue process is the /i-transform of the process 
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consisting of m independent BESQ(2(<5 — m + 1)) conditioned never to collide, as stated 
by Konig and O'Connell [15]. Another proof is given by Pcche ([21], page 68). Here, 
we will deduce the expression for q t (x,y) from p t (x,y) following Muirhead [20], namely, 
by projection on the unitary group: we will use the Weyl integration formula, then give 
a detcrminantal representation of hypcrgcometric functions of two matrix arguments. 
First, we state the Weyl integration formula Faraut [7] in the complex case: for any Borel 
function /, 



f(A)dA = C m / f(uau*)a{du)(V(a)ydai---da m , 

H m JU(m) JWL™ 

where C m = ir m ^ m -^ /T m (m), U(m) is the unitary group, a is the normalized Haar 
measure on U(m), a = diag(dj) and A = uau* . Hence, the semi-group of the eigenvalue 
process is given by James [12]: 

Qt(x,y) = C m (V(y) 2 ) / p t (x,uyu*)a(du) 

JU(m) 

„— (xi+Vi)/2t I TT„. \ / ip I t-xuyu 



T '\. V \V) ) TT „-(xi+vj)/2t I TT „. I ip , ... X 



(2t) m ( m +^T m (m)T m (m + v) .1^ 

where y — diag(j/j), a; is a positive definite matrix whose eigenvalues are x\, . . . , x m , qFi 
in the third line is a hypergeometric function with two matrix arguments Gross and 
Richards [10] and S = m + v. v > —1. Next, we need a lemma. 

Lemma 3. Let B,C € H m and let (6j), (cj) be their respective eigenvalues. Then 

p F q ((m + Hi)i<i< p , (m + <f>j)i<j< g ;B, C) 
= 7r ™(m-i)/2( P -g-i) rm ( m ) 

' (r(w + i)r 



n 



^1 T m (m + fii) 

T m (m + (j>j) det{p.<? q {{ni + l)i<j< p , (1 + <i>i)i<j<qMcf)i,f) 

ir 

3= 



n 



x ii (r(^ + 1))™ ^(s)y(c) 



/or aZZ fa, 4>j > — 1, 1 < i < p, 1 <j < q. 
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Proof. Recall that the hypergeometric function of two matrix arguments is given by the 
scries 

Miath^, (e 3 -)i<i< 9 , B, G) - n ? =1 (e,) T C T {I)k\ ' 

It is well known that 

Ct(-B) = 7— r L s r (&i, . . . ,6 m ), 
(m) T 

where s T is the Schur function and d T = s T (l, . . . , 1) is the representation trace or degree 
(Gross and Richards [10] and Faraut [7]). Substituting in the series gives 



p F q ((m + Hi)i<i< p , (m + 4>j) 



r> r \_'ST'Sr IlLl( m + ^)r S T {B)S T {C) 

i< j < q ,B,C)-^ UU(m + 4> . )T (m) T ■ 



We now write 

, , _ -^r T(m + m + k r -r+l) t r r(/i. t + 1 + k r + 5 r ) 

(m + (M) T 11 r( +m „ r + 1) -11 r( + m _ r + i) 

r— 1 v ' / r— 1 w 7 



* r ro (m + Mi ) M ( )fcr+5r ' 



m — r. 



Doing the same for each (m + <^j) T and for (m) T , we can see that 
p F q ((m + (H)i<i< p , (m + 4>j)x<j< q \ B, C) 



V^V^TTf nLi(^ + l)fc r +3,. \ s T (B)a T (C) 

^ ^ M v ni=i fa + / nr=i aw 



where 



/ m(m — 1) \ . „ . 

13= f-^ -J(p-9-1)- 

To obtain the desired result, we use the Hua formula Faraut [7]: 
Lemma 4. Given an entire function f{z) — X^feLo e k% t we have 

det(f{biCj))i d _ x^sT'l TT I s T (B)s r (C) 



ee n 



V(B)V(C) t^\r=\ ekr+&r ) s ^™) 



The Laguerre process and generalized Hartman-Watson law 
Thus, we obtain 



567 



p F q ((m + Hi)i<i< p , (m + 4>j)i<j< q ;B, C) 

= %m ( m -i ) ( P - q -i)/2 Tmim) Y[ r( Mi + 



LA T m (m + m) 

rWY" 00 nLi^' + 1 ) fc {hc P ) h \ 



11 i% + i) ^(B)V(CQ 



__.m(m-l)(p- 9 -l)/2 r / m \fr + -pr r m (m + 0j) 

m( ^^(m + ^/j r^- + i) 

det( p ^g((^ + l)i<i< p , (1 + </>j)i<j< q ; bic f )ij) 

x v(bmc) ' □ 

For p = and q > 1 , 

F q ((m + (f>j)i<j< q ]B, C) 

_„.-m(m-X)(g+X)/2 r / \ A £^(m+_^) det( J^((l + (j>j)\<j< q \ kc f )ij 
and similarly. 

pfRn _ r m (m) detfe 6 "*),,/ 
01 ' J TrmCm-^/a y(S)F(C) 

which can be viewed as a Harish-Chandra formula for the Itzykson-Zuber integral Collins 
[5]. We now proceed to the end of the proof. Taking p = 0, q=l, B = x/4t 2 , C = y, we 
obtain 

„( x \ (4t 2 ) m ( m - 1 )/ a r fB (m + i/)r wl (m)det(o^i((i/ + l);xitfi/4* 2 )) 



4i 2 ' W y ^m(m-X) ( r („ + !))m V(x)V(y) 

The expression for q t (x,y) follows from a simple computation and from the fact that 
oJ^i^ + l);^/^ 2 )^- _ / 2t V f Vx-y- \ 

i>+i) Kv^yJJ \ * )' 



□ 



Proposition 5. The measure defined by p(dx) = (det(x)) s m Ax on is invariant 
under the semi- group, that is, pPt = p. 
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Proof. Denote by P t the semi-group of the Laguerre process L(5,m,x) for 5 > m— 1. 
Then we have to show that 

P t f(x)p(dx) = [ f(y) P (dy), feC {H+). 
x>0 Jy>0 

This follows by a similar computation and the same arguments as in the proof of Corollary 
3. □ 

Remark 10. For Wishart processes, it is easy to see that 

p(dx) := (det(x)) s ^ m+1 ^ 2 l {x>0} dx 
is invariant under the semi-group. 

5. Girsanov formula and absolute-continuity relations 

The index v > — 1 of a L(S, m, x) is dehncd by v — S — m. In this section, we will proceed 
along the same lines as Donati-Martin et al. [6] to derive absolute-continuity relations 
between different indices. 

5.1. Positive indices 

Take a matrix-valued Hermitian predictable process H. Let Q s x be the probability law of 
L[5, m, x) for S > m — 1 and x G H^. Define 



-/() 



tv(H s dB s +H s dB s ) 



$ t =cxp^L t -i^ tr(ff s 2 )d S ^, 



where B is a complex Brownian matrix under Q x . We can easily see that the process f3 
defined by (3t = B t — L H s ds is a Brownian matrix under the probability 

Furthermore, (Xt)t>o is a solution of 

dX t = v^d/3 t + dp* ^X t + {\fx~ t H t + H t ^X t + 25I m ) dt. (7) 
For H t = v\[Xt 1 , (7) becomes 

dX t = ^X~ t dp t + d$Vx~ t + 2(5 + v)I m dt, 
so that (X t )t>o is L(S + u, m, x) under P x . Thus, we have proved the following theorem. 
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Theorem 6. For S > m — 1, 



v 
2 Jo 



fY, 1 dB< 



,2 ft 



2 Jo 



tr(X s _1 )ds 



Proposition 6. 



(8) 



(9) 



Proof. We know that V u (det(u)) = det(u)u _1 , hence, V tl (log(dct(w))) = vT 1 . Then, 
using the fact that, for S = m, (log(det(X t ))) is a local martingale, we obtain from the 
Ito formula that 

log(det(X t )) = log(det(X )) + / \x{Xr\VX,AB B + ^B k a ypQ ) ) 

Jo 

= log(det(X )) 



^XZ 1 dB, 



From (9), it follows that: 
Corollary 5. 



X t =y 



det(y)-» /2 p™ + »(x,y) _ T m (m) ,^,_^ v/i F 1 (m + v, z) _ I v (z) 



-(det(z))- 



oFi{m,z) i ( z )' 



det(x) PT( x iy) T m (m + v) 
where z = xy/At 2 . 

We now state the following asymptotic result: 
Corollary 6. Let X be a Laguerre process L(m,m,x). Then, as t — > oo 



tr(^)- 1 ds^T^P), 



(mlogt) 2 Jo 

where T\ is the first hitting time of 1 by a standard Brownian motion (3. 

Proof. From (9), we deduce that 
2v 2 



Q KA-(^l tT(X7l)ds 



X t =ty 



□ 
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2\ 



T m (m + 2v/m\ogt) oF^m^xy/it 2 ) 

Noting that (£™ 1 )- I V mlo g* — e _ly J and since both hypergeometric functions converge to 1 
as t — > oo, wc obtain 



(mlogi) 2 J 



J tv(X- 1 )ds\X t = ty 



Then, since 



p -tr(x)/2« 

lim t m p?(x, 2y) = lim e~ tr(a) fi m, ^ 

t— »oo t^oo L rn (m) 



-tr(y) 



2t / r m (m) : 



we obtain 
Q™fexp 



2v 2 



(mlogt) 2 



jf tx(Xj 1 )d^\x t = y^pT{x,y)Ay 



x exp 



(mlogi) 2 

2v 2 
(mlogt) 2 



J^triX: 1 )^ 



X t =ty)t m2 pT(x,ty)<iy t ^e 



by the dominated convergence theorem. 



□ 



5.2. Negative indices 

Take < a < dct(a;). The same computation as in Section 5.1 with H t = —v\fX~t , 
< v < 1 , shows that 

v / 2 / ,,2 rthT, 



?m — v\ 
X \^tAT a 



where T a := inf{t, det(X f ) = a}. Letting a — > and using the fact that To = oo almost 
surely under Q™ , we obtain 

_ / det(x) ' 
" \det(X t ), 



( det(x) ^ r\m+v\ 

Ix \&f 
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Proposition 7. For all t> and < u < 1, 

Proof. From the absolute-continuity relation above, we deduce that 



On the other hand, using the expression for the semi-group, 

<3S(det(*,n = ptr ^y ,Fi (-»;«;-£ 

. (2(r .E^) exp (_ tr (£)) lfl ( J + a;j; £ 

by the Kummer relation (cf. Theorem 7.4.3 in Muirhead [20]). Taking s = —v, we are 
done. □ 



6. Generalized Hartman— Watson law 

Henceforth, we will write & to denote one-dimensional hypergeometric functions. We 
define the generalized Hartman-Watson law as the law of 



/ tr^-^ds under Q™{-\X t =y). 
Jo 



Its Laplace transform is given by 

Q™(ex P (^ ttriX-^s) X t =y)=/f^ d et( Z )^^^, (10) 
\ \ 2 J / / T m (m + v) Q Fi(m,z) 

z = xy/4t 2 . Recall that for m = 1, this is the well-known Hartman-Watson law whose 
density was computed by Yor [23]. Here, we will investigate the case m = 2. The Gross- 
Richards formula is given, for p — and q = l, by 

Fi(m + v,z) = —— , 

V(z) 

where (zi) denote the eigenvalues of z and V(z) = Y\i<j( z i ~~ z j) i s the Vandermonde 
determinant. Noting that T m {m + v) = Iljli T( m + v ~ J : + 1), then 

= det(4 TO - j)/2 W^(2yi-)) 
det(4 m - i)/2 / m _ i (2V5D) ' 
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Without loss of generality, we will take t = 1. 

Proposition 8. For m = 2, let X\ > A 2 be the eigenvalues of ^Jxy. Then the density of 
the generalized Hartman-Watson law is given by 



/(«) 



p7tV27Tt; 3 

J °° zsinh(pVl - z 2 )c- 2 ^ T ^ zcosh ye^ 2 (y 2 - 7t2 '/ 1, (sinh 2 /)sin(47t 2 //v) dzdy 



Jo Jo ~ cosh j^^^ ) j Q (2^Ai A 2 mx) dudz 
/or f > 0, where p = Ai — A 2 . Furthermore, if Ai = A 2 := A, i/ien 

4A«e 2 " 2 /" / °° g (;/)e-^ 2 /"(sinhy) sin(47ty/t;) dy 

W 7t 2 VW i^ a (5;i;2;A 2 ) 

w/iere 

1 7T 7 2 (2Acosh?/) + L 2 (2Acosh?/) 

g{y) = g + o ^ — u ' 

3 2 2 A cosh j/ 

and L 2 is £/ie Struve function Gradshteyn and Ryzhik [9j. 

Proof. For m = 2, (10) becomes: 

, , = XJu+xJMIJM - A 2 /„+i(A 2 )/„(Ai) 
1 j A 1 / 1 (Ai)/o(A2)-A 2 /i(A 2 )/o(Ai) ' 

so, using the integral representations below (Brychkov et al. [3], page 46), 



x{al v j r \{ax)l v {bx) — bl u j r i{bx)l u {ax)) = {a —b) ul v {au)l v (bu)du 

Jo 

with x = 1, a = Ai, b = A2 and (Gradshteyn and Ryzhik [9], page 734) 



7t /a 



(I (>/^+fc))/,(f (>/^-6)) =/ acos y a r^ ) ^Md,, 

where a > 0, 9?(i/) > — 1. With a = 1, 6 = (Ai — A 2 )w :=pu and c = 2^/Ai A 2 u, the numera- 
tor of (10) is then equal to: 



-(A 1 -A 2 )/ / J2v(2VAiA 2 ux)duda;. 

71 Jo Jo VI -x 2 



Taking 1/ = 0, the denominator is then equal to: 

»i /.l 



-(A?-Aa) / / UC ° Sh{ ?;^h (2VMMu X )dudx. 
n Jo Jo Vl-x 2 
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Jo Jo U cosn (P M Vl — X 2 ) I y/l — X 2 I 2 u (2y / \iX 2 ux) du dx 



(10) 



Jo Jo wcosh(pu\/l — x 2 )/\/l — x 2 lofi^/ \1\2ux) dudx 
Now we only have to use the integral representation of l 2v (Yor [23] ) : 

2m J c 

2i7t7 c ; (2^3)V2 

where C is the contour indicated in Yor [23]; hence, the density function is given by 

^ ^_ 1 Jo 1 Jo Sc UUJ cosh (P u Vl - x 2 )/ Vl - xV^ 1 ^"* co^e- 2 " 2 ^ du da; du; ^ 
i-7r\/27Tt> 3 J 1 Jq 1 u cosh(pu\/l — a: 2 ) / \Jl — x 2 I (2VAiA 2 ux) du dx 

We can simplify this expression by integrating over C to see that the numerator is equal 
to (Yor [23]) 

\/AiA 2 f f 1 f 1 f°° ^ COSh(jWl ~ X 2 ) e _2VAlA?„xco S hfl 

7tV27Tu 3 Jo Jo Jo Vl — ^ 2 

x e- ( - 2te2 - 7t2))/u (sinhu)sin^^^ dudxdy. 

Setting z = ux, the numerator is written 

VAiA 2 u f 1 f u f°° ^ucosh^Vu 2 - z 2 ) ^_ 2VX ^ zcoshy 



7t\/ 27TU 3 7o Vo Vo Vu 2 — Z 2 



x c - 2 (« 2 -' r2 )/ t '(sinh 2 /)sin^^^ dudzdy, 
which wc can integrate with respect to u to obtain 

/' /"zsinh(pyT^)e- 2 ^ zcos ^e- 2 fe 2 -- 2 )A'(sinh y )sin^ N ld^du. 

p7TV27tW 3 JO JO V w / 

We now prove the second part. In this case p = and we have to evaluate 

Awe 2 " 2 /" Jo Jo J °° u 2 x/Vl - x 2 c- 2Au;rcosl ^c- 2 « 2 / t '(sinh 2 /) sin(47ty/u) dudxdy 
7TV2™ 3 £ f* ul (2Xux) / VT — a; 2 dit da; 
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Setting z = ux, the numerator reads 
Xve 2n2 / V f 1 f°° 



Integration with respect to z yields (Gradshteyn and Ryzhik [9], page 369): 



nV 2nv 3 

For the denominator, we use the fact that 



g(y)e- 2y2/ ^smhy)sm(^pj dy. 



which yields: 



— (zh(z)) = zl (z), 



• x ul (2\ux) = f 1 h(2Xx) ^ 



o Jo 



\/l — x 2 Jo 2Xx\/l — x 



Then, the formula 

x a - l {a 2 - x A f-'I v {cx) dx = 2~ v - l a 



T{p+(a + v)/2)T{y + \) 



. a + v a + u a c 

xA| + — ;^ + l; — ) 



taken with a = 0, a = 1,0= 1/2, c= 2A, v = 1 gives 

1 h(2Xx) , 7t (I „, 2 
dx = - 1 & 2 [ -;1;2;A 2 



lo 2XxVT^x 2 4 V2 

We can proceed differently. Let Ai = X 2 + h. Then 

((A 2 + h)I u+1 (X 2 + h)I v {X 2 ) - X 2 I V+X {X 2 )I V (X 2 + h))/h 



(10) 



((A a + h)h(X 2 + h)I (X 2 ) - A 2 /i(A 2 )/o(A 2 + h))Jh 



Next, we let h — > 0. As usual, we first compute the numerator and then take v = 0. To 
do this, we shall evaluate 

, ,. (X 2 + h)I v+1 {X 2 + h)-X 2 I u+ i(X 2 ) 
A = nm ; , 

h^Q h 

D .. I v (X 2 + h)-I v {X 2 ) 

B = Jim , 

h^O h 

which are equal respectively to ^(x/^+i(x)) and taken for x = X = X\ = X 2 . 

Using the differentiation formula -^{x v I v {x)) — x" I v -\{x) (Lcbedcv [18], page 110), we 
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thus 



■^-(xl v +i(x)) = -vl y+ i{x) + xl u (x), -^-(I v (x)) = --I v (x) +I„-i(x), 
ax ax x 



N = I v {X){-vI v+1 (X) + AI„(A)) - XI v +i(X) ( ™/„(A) + Iu-i(\) 



A 



^(A)-/ v+ i(A)J tf _ 1 (A) 



ig(A)-Ji(A)J_i(A) ' 
Using the integral representation (Gradshteyn and Ryzhik [9], page 757), 

2 f n/2 

I„{z)I v {z) = - / co^{ l i-v)6)I ll+l/ {2zcas6)de, »(/i + v)>-l 
71 Jo 

the numerator is written as 



N 



2 r' 2 

- / (1 -cos20)i" 2l ,(2Acos0)d0 



7t/2 

(sin 2 <9)/ 2l ,(2Acos0) d9 



\J\-r 2 I 2v (2\r) dr. 
Thus, using (11), the denominator is equal to 

D = ^J Vl^WXr) dr = ^ ^ Q; 1; 2; A 2 ^ . 
Finally, the integral representation of I v gives 

Awe 2 " 2 /" / °°e- 2 ^/ tl sinh(?/)sin(47T?;/w) J* rVl - r 2 e- 2Arcoshy drdw 



f(u) = 



nV2nu 3 Jg 1 VI -r 2 / (2Ar) dr 

Aue 2 " 2 /" J °° g(y)e~ 2y ~/ u sinh(y) sin(4m//it) du 

ttV2to 3 VT^7 2 / (2Ar)dr ' □ 
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Recall that, for < v < 1, 



«"-"< T » > <> = F^Sry *' (5) " m + - 5 



Proposition 9. Lei m = 2 and Ai > A2 6e ifte eigenvalues of x. The density of Sq := 
l/(2To) under Q™~ v is given by 

_ (\ 1 \ 2 yu 2 »- 2 C-^+ X ^ U ^(2, v+l, Am) - ifip, v+1, \ 2 u) 
HU) ~ T(v + l)T(v) (X1-X2) 

Corollary 7. If Xi = X 2 := A, the density is written 

2X 2l, u 2 "~ 1 c~ Xu 

/(u)= r(, + 2)i» + 

Proof. Recall first that when m = 1, So = Jv/x, where 7„ is a gamma variable with 
density r" _1 e _r dr. With the help of the Gross-Richards formula, it follows that, for 
m = 2, 

Q™-"(S <u) = (Ai yj =r u a "(Ai ^1(1/, 1/ + 2, -Am) - 1,1/ + 1, -A 2 u) 

(Ai — A2JI 2^ + 

- A 2 i^i(z/, v + 2, -A 2 u) 1^1 (1/ - 1,1/ + 1, -Aim)), 

where Sq := l/(2Tb). This is a C°° function in u. Hence, we will compute its derivative 
to obtain the density. Recall that 

^-x&i (a, b, z) = ^^(0 + 1,6+ l,z), 
dz b 



thus 



/(«) = ""(So < u ) = Ai, A 2 )w 2!/ - 1 (A - S), 



where 

(AiA 2 



if(i/, Ai,A 2 



r 2 (i/ + 2)(Ai-A a )' 
A = 2r/(Ai 1 ^1(1/, i/ + 2, -Aiu) 1^1 (1/ - 1, 1/ + 1, -A 2 u) 

-A 2 i^i(i/,i/ + 2,-A2«)i^i(i/-l,i/ + l,-Aiu)) ) 
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B = -^—(\\ui^i{v + + ^,-\ 1 u) - + 1,-X 2 u) 

- X 2 2 u 1^1(1/ + 1, v + 3, -X 2 u) 1^1 (1/ - 1, v + 1, -Xm)). 

Then, we use the contiguous relation 

b 1^1 (a, b, z) -bx^xia - 1,6, z) = Zi^i(a, b + 1, z) 

to see that 

Aim x &i{v + 1, 2/ + 3, -Xm) = {v + 2){ x &i{v, v + 2, -A x it) - 1 J*"i (1/ + 1, 1/ + 2, -Aim)), 
A 2 miJ ? iO + 1,i/ + 3,-A 2 m) = (i/ + 2)( 1 J?i(i/,i/ + 2,-A 2 u) - 1^1(1/ + l,i> + 2, -A 2 w)) 
imply that 

/(m) = Ki {v, A 1; A 2 )m 2 ^ 1 (C + D — E — F) 

where 

/fl(,/ ' Al ' A2) = r 2 (, + 2)(A 1 -A 2 )' 

C = Ai iJ^O, v + 2, -Aim) 1^(1/ - 1, v + 1, -A 2 m), 
D = Ai 1.^,(1/ + 1,1/ + 2, -Aim) - 1,1/ + 1, -A 2 m), 

E = A 2 i«^i(i/,i/ + 2,-A 2 u) i^i(i/ - 1,1/ + 1, -Aim), 
i^ 1 = A 2 1^1(1/ + 1,1/ + 2, -A 2 m) 1^1 (y - 1,1/ + 1, - Aim). 

Applying the above contiguous relation again yields: 

AiUiJ r i(i/ + l,2/ + 2,-A 1 u) = (1/ + l)(iJ*"i(i/,z/ + l,-Aiu) - i&i{v + l,i/ + l,-Aiu)), 
A 2 u !&i(v + l,u + 2, -A 2 it) = (i/ + l)(i^i(i/,z/+l,-A 2 u) - l,f + 1,-A 2 m)), 

A 2 k 1 J ? iO,2/ + 2,-A 2 m) = (i/+l)(i^i(i/- 1,i/+1,-A 2 m)- 1,-A 2 m)), 
AiMi^i(t/,i/ + 2,-AiM) = (i/ + 1)(i^i(i/- l,i/+ 1,-Aim) - i^i(i/,i/+ 1,-Aim)). 

Substituting in the expression for /, we obtain 

f{u) = K 2 {u,X 1 ,X 2 )u 2v - 2 {G-H), 

where 

t/( i / + l)(A 1 A 2 )^ 



K 2 (z/,Ai,A 2 ) 



r 2 (z/ + 2)(Ai-A 2 )' 
G = \&\{v + l,v+l,-X 2 u) \@\{y- l,u + 1,-Aiu), 
= 1^1 (u + l,v+ l,-X x u)x&i{v- l,u + 1,-A 2 u). 
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Eventually, writing 

T 2 {v + 2) = I> + 2)I> + 1) = v[y + 1)I> + 1)I», 
i^i (a, a, z) = e~ 2 i^i(a,6, -z) = e~ z \^\(b - a, 6, z), 

we obtain 

(AiA 2 )^ 2iy - 2 c-( Al + A2 ) u !^i(2, i/ + 1, Am) - i^i (2, v + 1, A 2 u) 



/(«) = 



r(^ + i)r(i/) Ai-a 2 



The case Ai = A2 is treated in the same way as before (for the Hartman- Watson law). 
In fact, writing Ai = A 2 + h and letting h — > 0, we see that the density is given by 

^ = I> + 1)I» dA^^^ + M.) 

= v( ^ i^i(3,i/ + 2, Au) 

2A 2,y ii 2,y - 1 c- A " 



i> + 2)l» 



-1^1 (z/- + 2, -Ait). 



□ 



8. Conclusion 

The Gross-Richards formula has been the main ingredient in this paper, since it enables 
us more explicitly to express the special functions with matrix argument. The case m = 3 
can be treated in the same way, but computation becomes too complicated. So, if we want 
to deal with the general case, it will be convenient to find a more explicit formula. Indeed, 
Schur functions can be expressed as polynomials in the elementary symmetric functions 
e r or as polynomials in the completely symmetric functions h r . More precisely, we have 

s\ = det(exi-i+j), l<i,j<n, 
sx = det(hy.-i+j), l<i,j<n, 

where A is a partition of length < n, and A' is the conjugate of A Macdonald [19]. So 
using these two identities, can we improve our results? 



Appendix: Special functions 
A.l. The hypergeometric series 



The multivariate hypergeometric functions were studied by Muirhead [20] in the real 
symmetric case, Chikuze [4] for the complex Hermitian case, and Faraut and Koranyi [8] 
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in a more general setting. For Hermitian matrix argument, they are denned by 

pFq((ai)i<i< p , (&j)i<j<«;^) =J2JL-St 7}TV~ CT b\ 

fe>OTii WT '" l(,?jT fc - 

where r = (fci, . . . , fc m ) is a partition of weight k and length m such that k\ > • •■ > fc m , 
(a) T is the generalized Pochammer symbol defined by 

c ^ TT r ( a + k i ~ 1 + 1 ) ru u ^ 

( a )r=M— ^ . , n , T={ki,...,k m ), 

and C T (X) is the zonal polynomial of X such that 

Several normalizations for this polynomial exist in the litterature, but we do not consider 
them. This polynomial is symmetric, homogeneous, of degree k in the eigenvalues of X, 
and is an cigenfunction of the differential operator 

- ft 2 m t 2 a 

jL —' Oxf t-* 1 Xi — Xh OXi 

t=l 1 1=1 l<kjti<m 

Furthermore, it is identified with the Schur function s T and C T {YX) = C r (VYX\^Y) 
for any Hermitian matrix Y. It is well known that, if p = q+l, then the hypergeometric 
series is convergent for < ||X|| < 1 (|| ■ || is the norm given by the spectral radius); if 
p< q, then it converges everywhere; otherwise it diverges. 

A. 2. The modified Bessel function 



The modified Bessel function Lebcdev [18] with index v £ K is given by the series 

oo s \ 2k+u 

It can be represented through standard hypergeometric functions o-Fi and iF\: 
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